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ABSTRACT 

Low-frequency observatories are currently being constructed with the goal of de- 
tecting redshifted 21cm emission from the epoch of reionization. These observatories 
will also be able to detect intensity fluctuations in the cumulative 21cm emission after 
reionization, from hydrogen in unresolved damped Lya absorbers (such as gas rich 
galaxies) down to a redshift z ~ 3.5. The inferred power spectrum of 21cm fluctua- 
tions at all redshifts will show acoustic oscillations, whose co-moving scale can be used 
as a standard ruler to infer the evolution of the equation of state for the dark energy. 
We find that the first generation of low-frequency experiments (such as MWA or LO- 
FAR) will be able to constrain the acoustic scale to within a few percent in a redshift 
window just prior to the end of the reionization era, provided that foregrounds can 
be removed over frequency band-passes of > 8MHz. This sensitivity to the acoustic 
scale is comparable to the best current measurements from galaxy redshift surveys, 
but at much higher redshifts. Future extensions of the first generation experiments 
(involving an order of magnitude increase in the antennae number of the MWA) could 
reach sensitivities below one percent in several redshift windows and could be used to 
study the dark energy in the unexplored redshift regime of 3.5 < 2 < 12. Moreover, 
new experiments with antennae designed to operate at higher frequencies would allow 
precision measurements (< 1%) of the acoustic peak to be made at more moderate 
redshifts (1.5 < z < 3.5), where they would be competitive with ambitious spectro- 
scopic galaxy surveys covering more than 1000 square degrees. Together with other 
data sets, observations of 21cm fluctuations will allow full coverage of the acoustic 
scale from the present time out to z ^ 12. 

Key words: cosmology: diffuse radiation, large scale structure, theory - galaxies: 
high redshift, inter-galactic medium 



1 INTRODUCTION 

Measurement of the fluctuations in the intensity of red- 
shifted 21cm emission from neutral hydrogen promises to be 
a powerful probe of the reionization era (Furlanetto et al. 
2006). The process of hydrogen reionization started with ion- 
ized (HII) regions around the first galaxies, which later grew 
to surround groups of galaxies. Reionization completed once 
these HII regions overlapped (defining the so-called overlap 
era) and fiUed-up most of the volume between galaxies. De- 
tection of the redshifted 21cm signal will not only probe the 
astrophysics of reionization, but also the matter power spec- 
trum during the epoch of reionization (McQuinn et al. 2006; 
Bowman et al. 2007). 

The conventional wisdom presumes that the 21cm sig- 



nal would disappear after the overlap epoch, because there 
is little neutral hydrogen left through most of intergalac- 
tic space. However, Wyithe & Loeb (2007) recently demon- 
strated that fluctuations in the 21cm emission would remain 
substantial over a range of epochs following the end of the 
overlap era owing to the significant fraction by mass of neu- 
tral hydrogen that is locked up in the dense pockets that 
form the damped Lya absorbers (DLAs) such as gas-rich 
galaxies. These systems trace the matter power spectrum on 
large scales. Hence observations of 21cm fluctuations could 
in principle be used as a cosmological probe both during the 
reionization era and in the post reionization IGM. 

The sky temperature, which provides the limiting factor 
in the system noise at the low frequencies relevant to 21cm 
studies, is proportional to (1 + z)^'^, and so is a factor of 
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~ 3.4[(1 + z)/5]^'® smaller at low redshifts than for obser- 
vations at 2: ~ 7. As a result, detectability of fluctuations in 
21cm emission may not decline substantially following the 
overlap epoch. When combined with the large fluctuations in 
redshifted 21cm emission during the reionization era, the de- 
tectability of a 21cm power spectrum after the end of reion- 
ization will allow the measurement of cosmological parame- 
ters over a wide range of redshifts. Much of this constraining 
power originates with redshift space distortions (McQuinn 
et al. 2006), which probe cosmology through the mapping 
between the vectors describing the wave number, and the 
emission frequency plus the visibility. The correct mapping 
produces an undistorted power spectrum. In this paper we 
consider baryonic acoustic oscillations (BAO). These pro- 
vide constraints on cosmology that are related to redshift 
space distortions, but which are particularly sensitive to the 
dark energy. 

The BAO scale provides a cosmic yardstick that can 
be used to measure the dependence of both the angular 
diameter distance and Hubble parameter on redshift. The 
wavelength of the BAO is related to the size of the sound 
horizon at recombination. Its value depends on the Hubble 
constant, and on the matter and baryon densities. However, 
it does not depend on the amount or nature of the dark en- 
ergy. Thus measurements of the angular diameter distance 
and Hubble parameter can in turn be used to constrain the 
possible evolution of the dark energy with cosmic time. This 
idea was originally proposed in relation to galaxy redshift 
surveys (Blake & Glazebrook 2003; Hu & Haiman 2003; 
Seo & Eisenstein 2003) and has since received significant 
theoretical attention (e.g. Glazebrook & Blake 2005; Seo & 
Eisenstein 2005; Seo & Eisenstein 2007; Angulo et al. 2007). 
Moreover, measurement of the BAO scale has been achieved 
within large surveys of galaxies at low redshift, illustrating 
its potential (Cole et al. 2005; Eisenstein et al. 2005). Galaxy 
redshift surveys are best suited to studies of the dark energy 
at relatively late times due to the difficulty of obtaining ac- 
curate redshifts for large numbers of high redshift galaxies. 
If the dark energy behaves like a cosmological constant, then 
its effect on the Hubble expansion is dominant only at z < 1 
and becomes negligible at 2; > 2. In this case studies of the 
BAO scale at low redshift would provide the most power- 
ful measurement. However, the origin of the dark energy is 
not understood, and so it is not known a priori which red- 
shift range should be studied in order to provide optimal 
constraints on possible theories for it. 

At high redshifts, baryonic acoustic oscillations in the 
21cm power spectrum should be detectable during the reion- 
ization era using future low frequency arrays (Bowman, 
Morales & Hewitt 2007). Following the end of the reioniza- 
tion era, the detection of a large number of individual galax- 
ies in redshifted 21cm emission could be used to trace the 
matter power spectrum using the future generation of radio 
telescopes (Abdalla & Rawlings 2005), in a manner entirely 
analogous to optical galaxy redshift surveys. On the other 
hand, 21cm observations using a low-frequency compact ra- 
dio array could also detect fluctuations in the total neutral 
hydrogen content within volumes of IGM dictated by the 
telescope beam and frequency band pass. Thus, in analogy 
to observations of the neutral IGM during the reionization 
era, one could construct the power spectrum of 21cm inten- 
sity fluctuations in the cumulative 21cm signal from all the 



unresolved pockets of neutral hydrogen in the IGM, regard- 
less of their mass. An important advantage of studying the 
matter power spectrum through detection of fluctuations in 
the total emission of unresolved sources, rather than in the 
space density of individual sources, is that the requirement 
of object detection is removed which, as we show, allows the 
power spectrum to be determined at much higher redshifts 
where individual sources may not be resolved in sufficient 
numbers. 

Our goal in this paper is to investigate the feasibility of 
using redshifted 21cm observations to make precise measure- 
ments of the scale of the BAO in the matter power spectrum. 
We show that observations of the acoustic scale using 21cm 
emission could be used to constrain the nature of the dark 
energy in the unexplored redshift range of 1.5 < z < 6, as 
well as during the reionization era, and so would be comple- 
mentary to galaxy redshift surveys. No other probes for pre- 
cision cosmology are currently being applied to this cosmic 
epoch (Corasaniti, Huterer, & Melchiorri 2007). The outline 
is as follows. In § [2] and [3] we will describe our approach 
for calculating the 21cm power spectrum, and its measure- 
ment uncertainties. The precision with which the scale of 
acoustic oscillations can be detected using upcoming 21cm 
observatories will be analyzed in § [J] with the correspond- 
ing constraints on the evolution of dark energy discussed in 
§ [5] Finally, we will summarize our main conclusions in § |6l 
Unless otherwise specified we adopt the set of cosmological 
parameters determined by WMAPS (Spergel et al. 2007) for 
a flat ACDM universe. 



2 THE POWER SPECTRUM OF 21CM 
FLUCTUATIONS 

A powerful statistical probe of the reionization era will be 
provided by the power spectrum of 21cm emission which 
is naturally accessible to interferometric observations such 
as those to be carried out by the Mileura Widefield Array 
(MWAQ) or the Low Frequency Array (LOFAeQ) . We may 
write the following expression for the power spectrum of 
21cm fluctuations 

P2l{k) ^b2l{z,kfP{k)D{zf, (1) 

where P{k) is the primordial power spectrum of the density 
field as a function of wave number k, extrapolated linearly 
to 2 = 0, and D{z) is the growth factor for linear pertur- 
bations. We model the power spectrum, including baryonic 
oscillations using the transfer function from Eisenstein & 
Hu (1998). The bias like quantity b2i{z, k), which has dimen- 
sions of temperature, may be determined from detailed nu- 
merical simulations or from an appropriate analytical model 
for bubble growth. The bias 621 will in general be a function 
of time and scale. 

We base our estimate for &21 on the modeling presented 
in Wyithe & Loeb (2007). At a specified redshift, our model 
yields the mass- averaged fraction of ionized regions, Qi, on 
various scales, R, as a function of overdensity, 5. We may 



1 See |http://www. haystack.mlt.edu/ast/arrays/mwa/index. html 
^ See |http: / /www.lofar.orgTl 
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Figure 1. Upper panel: The bias (f'21) computed from our semi- 
analytic model at 3 different scales (dark solid, dashed and dotted 
lines). For comparison the figure also shows three curves for the 
bias corresponding to the 21cm fluctuations given by spatially 
uniform mean neutral fractions, with fluctuations due to the den- 
sity field. In the three cases the neutral fractions are equal to 
1 (fully neutral, dashed grey line), 0.03 (the neutral fraction in 
DLAs at 2 < 5, dotted grey line), and the mean neutral frac- 
tion from our model (solid grey line). Lower panel: The bias (621) 
computed from our numerical model at 4 different wave numbers 
k (dark solid, dotted, short-dashed and long-dashed lines). The 
semi-analytic estimate at R = lOOMpc is also shown for compar- 
ison. 



then calculate the corresponding 21cni brightness tempera- 
ture contrast 



T{S, R) = 22mK 



1 + z 
7.5 



1/2 



[l-Qi{R,S)] ( 1+ 3<5) , (2) 



where the pre-factor of 4/3 on the overdensity refers to the 
spherically averaged enhancement of the brightness tem- 
perature due to pecuUar velocities in overdense regions 
(Bharadwaj & Ah 2005; Barkana & Loeb 2005a). Given 
the distribution of 5 from the primordial power spectrum 
of density fluctuations, we may find the probability dis- 
tribution dP/dT of brightness temperature T in redshifted 
21cm intensity maps. The second moment of this distribu- 
tion {(r — (T))^) corresponds to the auto-correlation func- 
tion of brightness temperature smoothed on a scale R. The 
effective bias can be estimated directly from our calculation 
of the auto-correlation function 

{{T-{T)f) 



621 (i?)' = 



(3) 



a{RY ' 

where (t{R) is the variance of the density fleld smoothed on 
a scale R. 

The upper panel of Figure [T] shows the bias (621) com- 
puted from our semi-analytic model. Curves are shown for 
model fluctuations computed at three different scales, each 
of which is larger than the typical bubble scale during the 
reionization era. For comparison, the flgure also shows three 
curves for the bias corresponding to the 21cm fluctuations 



given by spatially uniform mean neutral fractions, with fluc- 
tuations due to the density field. In the three cases, the neu- 
tral fractions are equal to 1 (fully neutral), 0.03 (the neutral 
fraction in DLAs at z < 5), and the mean neutral fraction 
from our model. We note that our analytic estimate of the 
bias is a weighted average of the bias on scales less than 
R, whereas the bias in equation ([1} should be computed 
at a particular scale. The bias computed using our analytic 
model is quite insensitive to scale for R > 30Mpc, indicating 
that on large scales the bias is independent of wave number 
k. 

Following the end of reionization (i.e. 2 < 6) there are 
no longer separate ionized bubbles as most of the IGM is 
ionized. Wyithe & Loeb (2007) have shown that the skew- 
ness of the 21cm intensity fluctuation distribution will be 
small during the post overlap epoch, which implies a lin- 
ear relation between the matter and 21cm power spectra. 
Moreover, since the mean free path of ionizing photons be- 
comes very large in the post overlap IGM, the ionization 
field is smooth. As a result the relation between the den- 
sity and 21cm emission can be reliably estimated using the 
semi-analytic model. 

However, during the reionization epoch the relation be- 
tween the power spectrum of 21cm fiuctuations and the 
underlying matter power spectrum is complex, and in the 
late stages is dominated by the formation of large ionized 
bubbles (Furlanetto et al. 2004). Thus prior to the end of 
reionization one must be careful about applying equation ((3]) 
to estimate the 21cm power spectrum in equation ((TJ. On 
very large scales such as those corresponding to the scale of 
baryonic acoustic oscillations, the fluctuations should aver- 
age over many bubbles so that the 21cm power spectrum 
is again expected to be linearly related to the matter power 
spectrum (McQuinn et al. 2006). Indeed, the largest discrete 
HII region that could be observed during the reionization era 
is much smaller than the BAO scale (Wyithe & Loeb 2004). 
Redshift space distortions can be ignored in measuring the 
acoustic oscillations scale during reionization, since the char- 
acteristic peculiar velocity is less than one percent of the 
Hubble velocity on this large scale. 

We have computed the matter and 21cm power spec- 
tra from a hybrid simulation (e.g. Messinger & Furlan- 
etto 2007) using both analytical and numerical techniques 
(GeU & Wyithe 2007). At each of redshifts z = 12, 2 = 8 
and z = 6.5 we computed the 3-dimensional ionization fleld 
within simulation boxes of side length 3000 co-moving Mpc. 
The simulations were computed with 256"^ resolution ele- 
ments in each case. For these simulations we used an input 
power spectrum which includes the BAO signal. At 2 = 12, 
2 = 8 and z — 6.5, the simulations have global neutral frac- 
tions of 98%, 48% and 11% respectively. The hybrid scheme 
of Geil & Wyithe (2007) is able to compute the neutral frac- 
tion within pixels that are larger than the typical size of HII 
regions. This feature allows us to compute power spectra at 
very large scales even early in the reionization era when no 
HII regions are resolved. 

The left hand panels of Figure[2]show the 21cm emission 
from 12Mpc slices through the simulation boxes at each red- 
shift. The higher redshift example (2 = 12) is early in the 
reionization era, and shows no HII regions forming at the 
resolution of the simulation (i.e. the IGM does not contain 
ionized bubbles with radii > 5 co- moving Mpc). The fluctu- 
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Figure 2. Examples of the 21cm power spectra during reionization. Left panels: Maps of the 21cm emission from slices through the 
numerical simulation boxes, each 3000 co-moving Mpc on a side with a thickness of 12 co-moving Mpc. In these maps yellow designates 
the absence of redshifted 21cm emission. Central-left panels: The corresponding matter power spectra multiplied by the growth factor 
squared (thin solid lines) as well as the 21cm (thick solid lines) power spectra computed from the simulation box. The input co-moving 
power spectrum P (also multiplied by the growth factor squared) is shown for comparison (short-dashed lines). Central-right panels: 
The baryonic wiggle component of the simulated 21cm power spectrum. The curves (thick dark lines) show the difference between the 
simulated 21cm power spectrum, and a reference no wiggle 21cm power spectrum computed from the theoretical no wiggle reference 
matter power spectrum multiplied by the square of the product between the bias and the growth factor [i.e. P21 — P^ef^%\D^\- For 
comparison, the red lines show the difference between the input matter and the no wiggle reference matter power spectra, multiplied by 
the bias and growth factor squared [i.e. (P — Pref )''2l^'^]' "^^^ grey band surrounding this curve shows the level of statistical scatter 
in realizations of the power spectrum due to the finite size of the simulation volume. Right panels: The scale dependent bias (621 )■ The 
upper, central and lower panels show results at 2 = 12, 2 = 8 and 2 = 6.5, which have global neutral fractions of 98%, 48% and 11% 
respectively in the model shown. 



ations in the 21cm emission are dominated by the density 
field at this time. The central redshift [z ~ 8) shows the 
IGM midway through the reionization process, and includes 
a few HII regions above the simulation resolution. The lower 
redshift example is just prior to overlap, when the IGM is 
dominated by large percolating HII regions, which are well 
resolved at the 5Mpc resolution of our simulation. 

The aim of this paper is to investigate the utility of 
21cm observations of the BAO scale as a probe of dark en- 
ergy at high redshift. It is therefore important to demon- 
strate that the percolation process does not wash out the 
signature of BAOs in the 21cm power spectrum towards the 
end of the reionization era. In the central left panels of Fig- 
ure [1] we show the corresponding matter [A^ = A;^P/(27r^)] 
and 21cm [A21 = k'' P21 / (2%^)] power spectra. The z = 6.5 



simulation exhibits a shoulder in the 21cm power spectrum 
at a scale corresponding to the characteristic bubble size. 
This shoulder is also seen in analytic models of bubble 
growth (Furlanetto, Zaldarriaga & Hernquist 2004). In the 
two higher redshift cases this shoulder is located at a scale 
below the simulation resolution. We may use these simula- 
tions to estimate the bias [621 (fc)] which is calculated from 
the square root of the ratio between the simulated 21cm 
power spectrum and the simulated matter power spectrum 
(using equation [3)1 . In the right hand panels of Figure [2] we 
show the resulting values of 621 (fc) as a function of scale. 

To illustrate the presence of BAOs in 21cm power spec- 
tra computed using these simulations, we have also cal- 
culated the difference between the 21cm power spectrum 
constructed from the simulation and a no wiggle reference 
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21cm power spectrum [i.e. P2i(fc) — f'2i,rcf (fc)]- This dif- 
ference is plotted in the central right panels of Figure [S] 
To construct P2i,ref (fc) we have multiplied a theoretical no 
wiggle reference matter power spectrum (Prof) by the bias 
computed from the same simulation box [i.e. P2i,rof(fc) = 
Prcf (fc)62i (fc)-D(z)^]. Note that since the boxes at each red- 
shift were generated using the same realization of the matter 
power spectrum, the noise is correlated between the result- 
ing 21cm power spectra. For comparison we also plot the 
difference between the theoretical matter power spectrum 
and the theoretical no wiggle reference matter power spec- 
trum, multiplied by the bias and growth factor squared [i.e. 
(P — Prof)&2ii^^, which is shown by the red line in Figure[2]. 
The grey band around this line illustrates the level of sta- 
tistical scatter in realizations of the power spectrum due 
to the finite size of the simulation volume (e.g. Peacock & 
West 1992). These simulations demonstrate that the 21cm 
power spectrum will exhibit BAOs throughout the reioniza- 
tion epoch including the percolation phase of HIl regions. 

Our simulations show that 621 will be constant on scales 
much larger than the characteristic size of HII regions, but 
scale dependent at larger values of k. The critical fc-scale 
below which 621 is nearly constant moves to smaller val- 
ues of k as reionization proceeds and the ionized bubbles 
grow in size. Early in the reionization process our results 
show that there is only a very weak dependence of 621 on 
scale since the fluctuations are driven by the density field. 
Indeed, our modeling shows the scale dependence of 621 to 
be weak on the scales of interest for BAO at all times until 
the end of the reionization era. At this time, just prior to 
the fuU overlap between HII regions, the simulations show a 
strong dependence of 621 on scale for k > 0.05 — 0.1 Mpc"^ 
(comparable to the BAO scale) due to the formation of large 
bubbles 0. Our simulations suggest that the scale dependent 
bias distorts the observed power spectrum of BAOs (see the 
panels in Figure [5] corresponding to z = 6.5). However this 
distortion, which causes the ratio of peak heights to increase 
relative to a spectrum with constant bias, is not expected 
to effect the extraction of the BAO scale. This is because in 
practice, the matter power spectrum could be fitted to the 
data using a scale dependent bias (Seo & Eisenstein 2005). 

The numerical values of &21 are plotted as a function of 
redshift in the lower panel of Figure [T] Also shown is the 
semi-analytic estimate of 621 on a scale of P = lOOMpc, 
which has been re-plotted for comparison. Our analytic and 
numerical models predict similar values of 621 and similar 
behavior with redshift. In particular, &21 has a value of a 
few tens of mK both at the very beginning and very end of 
reionization, despite the very different values of the global 
neutral fraction. In addition, on scales much greater than 
the typical bubble size there is a local minimum in the value 
of 621 midway through reionization, corresponding to the 
shift from fluctuations in 21cm emission being dominated 
by fluctuations in the density to fluctuations in the ioniza- 
tion field. However, the semi-analytic and numerical models 
do not agree in detail. The semi-analytic model does not 



Note that the bias may never be observed to be scale dependent 
at values of k smaller than the BAO scale since the evolution of 
the power spectrum becomes more rapid than the light crossing 
time of the BAO scale (Wyithe & Loeb 2004). 



include a Poisson component of fluctuation due to the finite 
number of bubbles in a region of radius R and so underes- 
timates the value of 621. On the other hand, our numerical 
scheme does not conserve photons [reflecting a limitation 
of semi-numerical models of this type (Messinger & Furlan- 
etto 2007)]. As a result, while our model predicts the topol- 
ogy of HII regions at a particular value of neutral fraction, it 
does not correctly predict the relation between the average 
neutral fraction and redshift. In the lower panel of Figure [T] 
this manifests itself as an offset in the redshift where the 
local minimum of the bias is predicted to occur. 

In light of the results described above, we have chosen 
to use the value for 621 computed at P = lOOMpc based 
on our semi-analytic model through the remainder of this 
paper. The semi-analytic model makes a conservatively low 
estimate of &21 at all redshifts and so will yield conservative 
estimates for the sensitivity of upcoming 21cm facilities to 
the 21cm power spectrum. 



3 MEASUREMENT UNCERTAINTIES IN THE 
POWER SPECTRUM OF 21CM 
FLUCTUATIONS 

Calculations of the sensitivity to the 21cm power spectrum 
for an interferometer have been presented by a number of 
authors. We follow the procedure outlined by McQuinn et 
al. (2006), drawing on results from Morales (2005) and Bow- 
man, Morales & Hewitt (2006) for the dependence of the 
array antenna density on radius, p{r). The uncertainty in a 
measurement of the power spectrum has two separate com- 
ponents. The first, due to the thermal noise of the instru- 
ment is 



AP2i,n(A:) 



Btint n{k±_) 



(4) 



where n{k±) is the density of baselines that observe the 
transverse component of the wave-vector (Morales 2005; 
Bowman, Morales & Hewitt 2006), T^ys ~ 250[(l-|-2)/7]^ '^K 
is the system temperatur^fl of the telescope when observing 
the 21cm line at redshift z, B is the band-pass over which 
the measurement of the power spectrum is made and tint 
is the integration time. The quantities D and AD are re- 
spectively the co-moving distance to the survey volume, and 
the co-moving depth of the survey volume (corresponding to 
the frequency band-pass within which the power spectrum 
is measured). The second component of uncertainty is due 
to sample variance within the finite volume of the survey, 
and equals 



(5) 



AP21,Sv(fc) =P2l(fc)^/— . 



The noise is evaluated within a fc-space volume ele- 
ment d^k. The total noise within a finite fc-space bin 
may then be obtained by integration over the volume 
within the bin. In both equations Q and ([5]) the quan- 
tity A^c = 27rfc^ sin6'dfcd6'[V/(27r)^] denotes the number of 
modes observed within a fc-space volume element d'^k = 



We assume Tsys to be dominated by the sky throughout this 
paper. 
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sm{6)dkd9. Note that in computing A'c we have as- 
sumed symmetry about the polar angle and expressed the 
wave vector k in components of its modulus k and angle 9 
relative to the line-of-sight. Because wave numbers can only 
be observed if their line-of-sight component fits within the 
observer's band-pass, we set A^c = if 2n/k cos{9) > AD. 
The number of modes observed depends on the volume of 
the survey, V — AD{\'^ /Auic), where Auic is the total 
physical surface area of an antenna (this point is discussed 
further below). 

The sensitivity to the 21cm power spectrum is depen- 
dent on both the sensitivity of the telescope to a particu- 
lar mode, and to the number of such modes in the survey. 
The former is set by the effective collecting area (Ac) of 
each antenna element (as well as the total number of an- 
tennae), while the latter is sensitive to the total physical 
area covered by each antenna (which we refer to as Atuc)- 
For a traditional interferometer consisting of a number of 
dishes in a phased array, these two areas are approximately 
equivalent Ac ~ Atuc, since the solid angle of the primary 
beam and the sensitivity are both proportional to the phys- 
ical collecting area of the dish. However in constructing the 
above formalism for the sensitivity to the power spectrum, 
we have explicitly allowed Ac 7^ ^tiic- This is because future 
interferometers being built to measure fluctuations from the 
epoch of reionization (like the MWA) will not be comprised 
of dishes, but rather a large number of 'tiles', each consisting 
of a phased array of A'dip dipoles distributed over an area 
^tiie- Since the size of the dipole will be much lower than 
A for observations of the 21cm line at 2: > 3.5, the effective 
collecting area of each tile in this regime is 
(Bowman et al. 2005) . Each tile forms an electronically steer- 
able primary beam, with solid angle nboam ~ A^/Atno. The 
MWA is designed to observe the 21cm line from the epoch 
of reionization, and so has {Ac ~ ^tiio) when observing at 
2: ~ 8, but Ac < ^tiic at lower redshiftijfl- In terms of mea- 
suring a power spectrum this reduces the efficiency of the 
MWA at higher frequencies. In this paper we consider the 
power spectrum at 1.5 ^ z < 6 as well as during the epoch 
of reionization ai z > 6. When showing results at z > 3.5 we 
use the specifications of the MWA since this observatory is 
already under construction and its design is not flexible. Ob- 
servations of neutral hydrogen at z < 3.5 are not accessible 
to the MWA, and so a new telescope would need to be con- 
structed for probing this epoch. Thus, at z < 3.5 we assume 
Ac ~ ^tiio, corresponding to a telescope with an optimal 
design (i.e. with dipoles spaced by A/2) for observations at 
these lower redshifts. 

In the case of a spherically averaged power spectrum, 
P2i(fc), each of the above noise components can be com- 
puted within fe-space volumes of d^k = 47rfe^Afe, where Afc 
is a finite bin of values in k. However, the power spectrum, 
f2i(fc±,fc||), can also be expressed in terms of the wave 
vector components that are parallel (fcy) and perpendicu- 
lar (fcx) to the line-of-sight, in which case the components 
AP2i,N and AP2i,sv can be computed within fc-space vol- 
umes of Sk — 27rfcx Afc^ Afc|| . 

The contamination by foregrounds provides an addi- 



^ The relation between Ac and Ax_i\c is often expressed in terms 
an aperture efficiency e = Acj At^i\c (e.g. Morales 2005). 



tional source of uncertainty in the estimate of the power 
spectrum. McQuinn et al. (2006) have shown that it should 
be possible to remove the power due to foregrounds to a 
level below the noise in the cosmological signal, provided 
that the region of frequency band-pass from which the power 
spectrum is estimated (B) is substantially smaller than the 
total band-pass available (Btotai)- Following the approxima- 
tion suggested in McQuinn et al. (2006), we combine the 
above components to yield the uncertainty in the estimate 
of the power spectrum. For the spherically averaged power 
spectrum we assume 



AP2l(fe) 



AP21,Sv(fc) + AP21,N(fc) 



t/S 



if fc > fc mill 

otherwise 



(6) 



while for noise in the power spectrum at k — (fcy , fej 
take 



AP2l(fc) 



AP2I Sv(fc) + AP21.N(fc) 
V^ficldB Stot/S 



if ^11 ^ ^min 

otherwise 



(7) 



where kmin = 2n/AD. In each of equations (O and (O 
the denominator represents the number of independent mea- 
surements made of the power spectrum within a band width 
B. The factor A'fleids is included because one independent 
measurement would be made per field imaged for time tint. 
As part of their analysis of foreground removal, McQuinn 
et al. (2006) have found that foreground power could be re- 
moved within a region B of the observed band-pass provided 
B is significantly smaller than Ptot . However in addition Mc- 
Quinn et al. (2006) have also found that foreground removal 
is not sensitive to the location of B within the total pro- 
cessed band pass. As a result, on scales where the power 
spectrum can be measured (i.e. at fc > fcmin) it can be de- 
termined within Btotai/P independent regions of the total 
processed band pass Ptotai- Hence while foreground removal 
will limit the scale of fluctuations that can be observed, fore- 
grounds should not affect the total sensitivity of the array 
to smaller scale modes. We note that, as with all studies of 
21cm fluctuations, the level to which it will be possible to 
remove foregrounds provides the largest uncertainty in our 
analysis. 



4 21CM OBSERVATIONS OF BARYONIC 
ACOUSTIC OSCILLATIONS 

Figure [3] shows results for the power spectrum of 21cm fluc- 
tuations, and sensitivity to the BAG scale at 2 = 3.5, 5, 
6.5 and 8. The spherically averaged model power spectra 
are marked on the left-hand panels as the thick dark lines 
(these are only valid on large scales as they do not capture 
the scale dependent bias that is the signature of the HII 
bubbles on smaller scales). We also show (heavy grey lines) 
the component of the power spectrum due to the baryonic 
acoustic oscillations. (Note that we show the absolute value 
of the full power spectrum minus a representative no wiggle 
power spectrum.) Estimates of the sample variance (dotted 
lines) and thermal noise (dashed lines) components of the 
uncertainty for detection by the Mileura Widefield Array 
(MWA) are plotted in each of the panels. The MWA, which 
is currently under construction will comprise a phased array 
of 500 tiles. Each tile will contain 16 cross-dipoles to yield an 
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Figure 3. The power spectrum of 21cm fluctuations and measurement of the baryonic acoustic oscillations. Results are shown at four 
redshifts, z = 3.5, 5, 6.5 and 8, and assume a low-frequency array with the specifications of the MWA, lOOOhr of integration on a 
single field and foreground removal within B = 8MHz of band-pass. Left panels: The power spectrum A^(fc) = A;^P2i(fc)/(27r^) of 21cm 
fluctuations (thick solid lines). The absolute value of the component of the 21cm power spectrum due to the baryonic acoustic oscillations 
is also plotted at fc < O.SMpc"^ (heavy grey lines). Also shown for comparison are estimates of the noise. In each case we plot the sample 
variance (dotted lines) and thermal noise (dashed lines) components of the uncertainty within fc-space bins of size Ak = k/10. The 
combined uncertainty including the minimum k cutoff due to foreground subtraction, is also shown as the thin solid lines. Central panels: 
The power spectra with the representative smooth power spectrum subtracted. The points with error bars show the accuracy attainable 
within a bin of width Ak/k = 0.1. The vertical dotted line is the wave number corresponding to the band-pass, below which the error 
bars are very large. Note that the vertical scale is different at each redshift and has been chosen to best illustrate the magnitude of the 
uncertainty relative to the amplitude of the BAOs. Right panels: The probability distributions for the recovered acoustic wave number 
k\ at each redshift. 



effective collecting area of = 16(AV4) (the area is capped 
for A > 2.1m). The tiles will be distributed over an area with 
diameter 1.5km. The physical area of a tile is Atue = 16m^. 
In this paper we consider 1000 tiles for phase II of the MWA. 
We model the antennae distribution as having p(r) oc 
with a maximum radius of 750m and a finite density core of 
radius 20m, and we assume a lOOOhr integration on a sin- 
gle field, a band-pass over which foregrounds are removed 



of B = 8MHz, and fc-space bins of widtljfl Afc = fc/lO. The 
total processed band pass for the MWA is Btotai ~ 32MHz. 
The combined uncertainty including the minimum k cutoff 
due to foreground subtraction is shown as the thin solid line. 

The central panels of Figure [3] show the power spec- 
tra at z = 3.5, 5, 6.5 and 8 with the representative smooth 



^ The signal-to-noise is increased in proportion to V Ak, and so 
will be substantially better per bin in measurements of the power 
spectrum at lower resolution in k. 
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power spectrum subtracted. The points with error bars show 
the accuracy attainable within a bin of width AA:/fc = 0.1. 
The vertical dotted line is the wave number corresponding to 
the band-pass, below which the error bars are very large. We 
have fitted the analytic approximatioiQ to the baryonic os- 
cillation component of the spherically averaged power spec- 



More recently a new technique has been proposed (Percival et 
al. 2007; Angulo et al. 2007) which replaces the analytic form of 
equation with a scheme that uses a reference power spectrum 
derived from the observed power spectrum, a full linear pertur- 
bation theory power spectrum, plus a damping scale to account 
for non-linear evolution. This method improves the fit to BAO 
power spectra computed in numerical simulations, and provides 
a more general approach. However we have chosen to employ the 
simpler approach of Blake & Glazebrook (2003) in this initial 
investigation. 



trum following Blake 
P(fc) 



Prcf(fc) 



= 1 -f Ak exp 



Glazebrook (2003) 
k 



0.07Mpc" 



27r ^ 



(8) 

This function has two parameters A and kA- The value of 
A is determined to high accuracy from observations of the 
CMB. For the purposes of this analysis we therefore assume 
that ^ is a known constant, and fit only for kx- We fit 
only to values of fc < 0.25Mpc~^. The accuracy to which 
kA can be measured determines the constraints that 21cm 
power spectra can place on the dark energy. The right-hand 
panels of Figure [3] show the probability distributions for the 
recovered /ca at each redshift considered. In a single field the 
MWA could detect the acoustic scale just prior to overlap, 
but could not make a precise measurement (less than a few 
percent) at any redshift. 

At values of fc ^ 10~^ Mpc"'^, the measurement of the 
power spectrum using the MWA will be limited by the ther- 
mal sensitivity of the array, and so the signal-to-noise achiev- 
able in this regime will be greatly enhanced by a subsequent 
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Figure 5. The relative accuracy achievable on the measurement of k\ as a function of redshift, assuming a prior probability on k\ of 
kA > 0.03Mpc-l. Results for the MWA (left panel) and the MWA5000 (right panel), with foreground removal in S = 6MHz (dotted 
lines), B = 8MHz (solid lines) and B = 12MHz (dashed lines) band-passes respectively. The upper and lower panels show results assuming 
1 and 3 fields, respectively (with each field observed for lOOOhr). In each panel we plot a vertical line at z = 3.5. To the right of this line 
we assume antennae with the specifications of the MWA. To the left of this line we assume an antennae number equal to the MWA or 
MWA5000, but tiles with Ae = Atue- 



generation of instruments with a larger collecting area. As an 
example, we consider a hypothetical follow-up instrument to 
the MWA which would comprise 5 times the total collecting 
area. We refer to this follow-up telescope as the MWA5000. 
The design philosophy for the MWA5000 would be similar to 
the MWA, and we therefore assume antennae distributed as 
p(r) oc with a diameter of 2km and a flat density core of 
radius 80m (see McQuinn et al. 2006). In Figure U we repeat 
our analysis of the power spectrum and baryonic acoustic os- 
cillations at z = 3.5, 5, 6.5 and 8 for measurements using the 
MWA5000. The panels show the same results as described 
in Figure [3l For the model overlapping at z = 6, we find that 
the BAO scale could be detected at at a range of redshifts 
and that very good measurements of /ca (~ 1%) could be 
made at z = 6.5 with the MWA5000 in a single field. 

Following the results of McQuinn et al. (2006) we do 
not make estimates for the SKAp in this paper. Current 
projections for the specifications of the SKA call for large 
antennae, with a small fraction of collecting area concen- 
trated in a core. This design limits the field of view, as well 
as the fraction of the telescope that can be used to measure 
the large scale modes which probe the BAOs. Thus, despite 
its increased collecting area, the SKA would be less powerful 
(with respect to measurement of the redshifted 21cm power 
spectrum) than the MWA5000, whose design would be opti- 
mized for the measurement of the 21cm power spectrum at 
high redshift. We note that at some redshifts the MWA5000 



* See www.skatelescope.org/ 
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would be cosmic variance limited on scales relevant to BAO 
studies. As a result if an SKA were built with a design based 
on the MWA5000, but with ten times the collecting area, no 
substantial gains could be made using observations of an in- 
dividual field. Of course a telescope with a larger collecting 
area could reach the limit of cosmic variance in a shorter 
integration, allowing more fields to be observed. 

To quantify the relative accuracy achievable on the mea- 
surement of kA we plot the variance of the recovered distri- 
bution divided by the best fit value in Figure[Sl In each panel 
of Figure [5] we plot a vertical line at z = 3.5. The MWA and 
MWA5000 could observe to the right of this line, and we 
assume antennae with the specifications of the MWA in this 
region. We note that there is non-zero probability for the 
recovered to lie at a harmonic of the true value (since 
the fitting function is quasi-periodic and the data is noisy). 
This may be seen in the distributions shown in Figures [3]- 
(4] For the results presented in Figure [5] we assume a prior 
probability on kA which is constant for kA > 0.03Mpc~^ 
but equal to zero otherwise. This effectively assumes that 
we know the accuracy of kA to ~ 20% apriori. The upper 
left hand panel of Figure [5] shows results for the MWA, with 
B = 6MHz, B = 8MHz and B = 12MHz band-passes re- 
spectively, and lOOOhr of integration for a single field. If 
foreground subtraction could be achieved, the larger band- 
passes would improve the accuracy significantly at low z 
by giving access to the peak centered on fc ~ 0.05Mpc~'^. 
The MWA will not make a precision measurement of the 
BAO scale using only one field. The upper right hand panel 
of Figure [5] shows the corresponding results for MWA5000, 
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again for a single field and lOOOhr of integration. In a sin- 
gle field the MWA5000 could make precise measurements 
(A/cA/fcA ~ 1%) over an extended interval prior to overlap, 
but not at higher redshifts. 

Since any one observing field can only be observed for a 
fraction of the time, measurement of the 21cm power spec- 
trum will be performed over several different fields. In addi- 
tion, some phased arrays will have the capability to observe 
using several primary beams at once. In the lower panels of 
Figure [5] we show results that assume an integration time 
of lOOOhr on each of 3 separate fields (note that this corre- 
sponds to ~ 3hr per day per field for 1 year). The noise on 
the power spectrum scales as the inverse square-root of the 
number of fields (McQuinn et al. 2006), which results in an 
improved precision on measurements of /ca over that achiev- 
able in a single field [see equations (|6} and ((7}]. The preci- 
sion achieved at z ~ 6.5 would be as low as Akx/kA ~ 3% 
using the MWA (comparable with the best current mea- 
surements from galaxy surveys), while the MWA5000 could 
reach A/cA/fcA ~ 0.5%. In addition the MWA5QQQ could pre- 
cisely (AfcA/AiA ^ 2%) measure the acoustic scale at z ~ 3.5. 

The antenna design of the MWA is optimized for the 
epoch of reionization measurements and only allows the 
21cm line to be observed at z > 3.5. However, in Figure [5] 
we show results down to a redshift of z = because a fu- 
ture instrument could be constructed with antennae that 
are sensitive to a different frequency range. At z < 3.5 (to 
the left of the vertical grey lines in Figure [5)l we assume 
Ac = Atiic and that the system temperature is dominated 
by the sky (an assumption implying that, unlike the MWA, 
such a telescope will need to have cooled receivers due to 
the lower sky temperature at shorter wavelengths). An in- 
strument constructed with the same number of antennae as 
the MWA but which operated at a higher frequency range 
could accurately measure the scale of acoustic oscillations at 
lower redshifts provided that foreground subtraction could 
be achieved over a sufficiently large band-pass. In construct- 
ing Figure [S] we have computed the 21cm power spectrum 
assuming a constant value of &21 at « < 2. Our modeling of 
621 does not include quasars and assumes local absorption of 
ionizing photons. It therefore becomes unreliable at 2 < 2. 
However the assumption of constant 621 is reasonable given 
that the observed density parameter of neutral gas does not 
vary significantly with redshift (Prochaska et al. 2005). At 
low redshifts the precision for measurements of ^a would be 
limited by the ability to remove foregrounds. The minimum 
value of k which is probed by the 21cm power spectrum be- 
comes larger as the observing frequency is increased, and 
at 2 < 1 has moved to a value beyond the scale of the 
BAO peaks. Nevertheless, at z ~ 1.5 observations from 3 
fields could yield precisions of ~ 2% and ~ 1% on mea- 
surements of feA using 1000 and 5000 antennae (Ac = Atiic) 
respectively, with even higher precision (~ 0.7% and ~ 0.5% 
respectively) attainable at z ~ 2.5. 



4.1 Sensitivity to the transverse and line-of-sight 
acoustic scales 

In the previous section we have computed the sensitivity 
of 21cm experiments to the angle-averaged value of fcA. 
However, observations of the 3-dimensional power spectrum 
of 21cm fluctuations provide constraints on both the ra- 



dial and transverse measures of this scale. In Figure [6] we 
present results for the sensitivity of 21cm observations to 
the line-of-sight and transverse BAO scale &t z = 3.5, 
z = 5, 2 = 6.5 and 2 = 8. In the left hand panels we 
show the signal-to-noise for observations of the full power 
spectrum P2i(fc±,fc||). In the central panels we show the 
corresponding signal-to-noise for observations of the differ- 
ence between the full power spectrum P2i(k±,k\\) and the 
no wiggle reference power spectrum P2i,rcf(fc±, fc||)- For this 
calculation we assume observation of a single field using 
the MWA5000, lOOOhr of integration, and foreground re- 
moval within B = 8MHz of band-pass. The signal-to-noise 
has been computed in bins of volume 27rfcx Afcx Afcy , where 
Afcx = fcx/10 and Afcy = fc||/10. We have used the analyti- 
cal approximation from Glazebrook & Blake (2005) 



P{k) 

Prof(fc) 



1 -f Ak exp 



27r 



0.07Mpc" 



«A,X 



+ 



fc|| 
^A, 



(9) 



where k^ = fcx + fcy, to estimate the corresponding con- 
straints on the line-of-sight and transverse acoustic scales 
(/caji and fcA,x)- On the right hand panels of Figure |6] we 
show contours of likelihood for the recovered values of feA,x 
and fcA,|| around the true input value for our standard cos- 
mology. Figure |6] shows that redshifted 21cm observations 
would be sensitive to both the transverse and line-of-sight 
components of the oscillation scale, with comparable uncer- 
tainty. 



4.2 Comparison with galaxy surveys 

Constraints on the BAO acoustic scale may be used to con- 
strain parameters in models of dark energy, with the abil- 
ity of a survey to discriminate among different models of 
dark energy governed by the accuracy achieved in measure- 
ments of the line-of-sight and transverse acoustic scales. Be- 
fore proceeding to discuss dark energy constraints, we there- 
fore first pause to compare the accuracy of the 21cm exper- 
iment with potential galaxy redshift surveys. Glazebrook 
& Blake (2005) present the simulated precision on mea- 
surements of AiA.x and kx,\\ from hypothetical galaxy red- 
shift surveys. Spectroscopic surveys of ~ 10^ galaxies within 
A2 = 0.5 redshift bins in the range 1 < 2 < 3.5 covering 
1000 square degrees would each measure the transverse and 
line-of-sight acoustic scales to accuracies of ~ 1% and ~ 2% 
respectively. On the other hand a photometric redshift sur- 
vey over 2000 square degrees between 2.5 < 2 < 3.5 would 
measure the transverse scale to ~ 1%, but would not con- 
strain the line-of-sight scale. Similar results were obtained 
at 1.5 < 2 < 2.5. 

The SKA could also be used to do a galaxy survey of 
sufficient size to measure the BAO scale (Abdalla & Rawl- 
ings 2004), so long as it were designed to have a sufficiently 
large field of view. However, the full sensitivity of the SKA 
would be required to push the galaxy survey beyond 2 ~ 1.5, 
limiting the studies of BAO using galaxies to relatively low 
redshifts when compared with optical spectroscopic surveys. 
By looking at fluctuations in the surface brightness of un- 
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Figure 6. Left panels: The signal-to-noise ratio (S/N) for observations of the power spectrum P2i{k_L,k\\)- Central panels: The signal- 
to-noisc ratio (S/N) for observations of the difference between the full power spectrum P2i(fcx,fc||) and the no wiggle reference power 
spectrum Prof,2l 1 )• In each case the signal-to-noise ratio has been computed in bins of volume 27rfcx Afcj^ Afcy , where Afcperp = 
k±/10 and Afcy = fcy/lO. Right panels: Contours of constant likelihood for the recovered fcA,± and fcA,|| around the true input value. 
Contours are shown at values of ~ Xmin ~ ^' 2.71 and 4, where Xmin value corresponding to the best fit parameter set. When 

projected onto individual parameter axes (fcA,x a-nd ''^A,||)i the extrema of these contours represent the 68%, 90% and 95% confidence 
intervals on values of the individual parameters. The four rows show results at z = 3.5, z = 5, z = 6.5 and z = 8. The observational 
parameters assume the layout and collecting area of the MWA5000 with lOOOhr of integration on a single field, a band-pass of 8MHz 
and redshifts of z = 3.5, 2 = 5, 2 = 6.5 and 2 = 8. 
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resolved 21cm emission, an array like the MWA could push 
measurement of the BAO to much higher redshift. 

The very large areas of sky that must be surveyed in 
order to measure the BAO scale using galaxy redshift sur- 
veys arise because very large volumes must be sampled in 
order to beat down the statistical noise on the large scale 
modes relevant to BAOs. For example, in units of the SDSS 
survey volume (Vsioan ~ 5.8 x lO^Mpc^), V/Vsioan ~ 3 and 
V/Vsioan ~ 1-8 are required to achieve 2% accuracy on kx at 
2 ~ 1 and z ~ 3 respectively (Blake & Glazebrook 2003). By 
comparison, in a single pointing within -Btotai = 32MHz of 
band pass the MWA surveys V/Vsioan ~ 2.4, V/Vaioan ~ 4.8, 

V/Vsloan ~ 7.6 and V/Vsloan ~ 11 at 2 = 1.5, z = 2.5, 

a = 3.5 and z = 6.5 (where we have assumed = Atwc) 
respectively. These volumes must be multiplied by A^floids in 
order to get the total volume from which the power spec- 
trum is to be constructed. Thus observations of the 21cm 
power spectrum will probe very large volumes of the IGM, 
comparable to the most ambitious galaxy redshift surveys. 

We find that observations using a low frequency array 
the size of the MWA, but with an appropriate frequency 
range would achieve measurements of the acoustic scale at 
z ~ 1.5-2.5 that are comparable in precision (~ 1%) to high 
redshift galaxy surveys using the next generation of opti- 
cal instruments. Moreover an instrument with the collecting 
area of the MWA5000 would extend this sensitivity out to 
z ~ 3.5. Observations of 21cm fluctuations would also allow 
the acoustic scale to be measured with comparable precision 
at much earlier cosmic epochs (z ~ 6), using the MWA5000. 
Moreover, 21cm observations will be comparably sensitive to 
the line-of-sight and transverse acoustic scales. In contrast, 
spectroscopic galaxy redshift surveys are more sensitive to 
the transverse scale, while photometric redshift surveys will 
not be sensitive at all to the radial BAO (Glazebrook & 
Blake 2005). 



5 CONSTRAINTS ON DARK ENERGY 

The measurement of the angular scale of acoustic oscillations 
(i.e. oscillations transverse to the line-of-sight) provides a 
measure of the angular diameter distance at a redshift z 
where the oscillations are observed. The angular diameter 
distance can also be computed at a redshift z for a chosen 
cosmology 



Da = {l + z)-^ 



dz' 

Hi^y 



(10) 



where H{z) is the Hubble parameter at redshift z, and we 
have set the speed of light equal to unity. Thus the mea- 
surement of -Da (2) probes the underlying geometry of the 
universe. On the other hand, measurement of the redshift 
scale of acoustic oscillations (i.e. oscillations along the line- 
of-sight) provides a direct measure of the Hubble parameter 
at a redshift z, which may be written 



H{z) = -Ho{l + z 



,3/2 



Sim. + SIq exp 3 



dz' 



(TTIO"'^^^'^ 



-,1/2 



In writing the Hubble parameter we have parameterized the 
equation of state as. 



Pi = Wi{z)pi 



(12) 



where the subscript i refers to either pressureless matter 
(with w = and a subscript i = m) or dark energy (with 
a subscript i — Q). We have also used the energy conser- 
vation for the dark energy pg = ~3H{z)[l + ■wq{z)]pq. 
The redshift dependence of the acoustic scale has been sug- 
gested as a powerful probe of the evolution of the dark en- 
ergy (Blake & Glazebrook 2003; Hu & Haiman 2003; Seo & 
Eisenstein 2003; Glazebrook & Blake 2005; Seo & Eisenstein 
2007; Angulo et al. 2007). If the size of the sound horizon 
is known, then the relative errors AD a/ Da in the angular 
diameter distance and AH/H in the Hubble parameter are 
related to the relative errors AA;A,x/fcA,± and AfcA,||/fcA,|| in 
the transverse and line-of-sight acoustic scales respectively 
(Figure O. 

To investigate the utility of the redshifted 21cm emis- 
sion as a probe of the dark energy equation of state we adopt 
a simple approach. In the remainder of this section we as- 
sume that the values of ilm, and the present-day value 
of the Hubble parameter Hq have been precisely determined 
apriori, so that we do not consider the uncertainty in the 
sound horizon. We also assume that the Universe is flat. In 
the upper panels of Figure [7] we plot the fractional change 
in the angular diameter distance and Hubble parameter re- 
spectively for several models with constant wq, relative to 
a model with a cosmological constant {wq = — 1). The ver- 
tical bars illustrate the type of precision achievable by the 
MWA (grey bar) and MWA5000 (dark bars) in a lOOOhr 
observation of three fields and a single fleld respectively, 
at the redshifts where these observations will be most ef- 
ficient. Based on the transverse measurement of the BAO 
scale, the MWA would be able to make a measurement of 
WQ with ~ 20% accuracy using a single measurement at 
z ~ 6.5, while the MWA5000 could make a measurement of 
Wq at the ~ 7% level from observations at each of z ~ 3.5 
and z ~ 6.5, yielding a combined constraint of better than 
5%. Better precision could be achieved using observations 
of three fields. Since a cosmological constant does not affect 
the dynamics at high redshift the line-of-sight BAO scale is 
left unaffected by this change in the equation of state. Thus 
a cosmological constant would introduce an asymmetry in 
the BAO scale that would be easily detected by high redshift 
observations. 

We can obtain a more quantitative relation between the 
uncertainties in feA and wq. If we expand around wq = — 1, 
then in the case of a constant wq model we can relate the 
uncertainty in wq to the observed uncertainty in the scale 
kA- We obtain 



AwQ — 



dwQ 



(11) 



d5,A ^^^^ 

where 5^1 = AkA,±/kA,± or 5^°/ = AfcA,||/fcA,|| for the 
transverse and line-of-sight components respectively. In the 
lower left panel of Figure[7]we plot the resulting uncertainty 
in Wq per unit fractional uncertainty in the transverse mea- 
surement of kA- This figure shows that a 10% uncertainty in 
Wq requires ^a to be measured in the transverse direction 
with a precision of ~ 1.5%. As mentioned above this will be 
achievable at each of redshifts 2 ~ 3.5 and 2 ~ 6 for the 
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Figure 7. Constraints on a constant wq. Upper panels: The relative cliange in angular diameter distance (left panel) and Hubble 
parameter (right panel) as a function of redshift for different models of the evolution in wq. The fractional changes are relative to a 
cosmological constant {wq = —1). The thick black and grey bars illustrate an accuracy of 1% and 3% on kj^ which would be achievable in 
some redshift ranges in a lOOOhr observation by the MWA5000 (1 field) and MWA (3 fields) respectively. Lower panels: The uncertainty 
in WQ per unit relative uncertainty in the angular diameter distance Da (left panel) and Hubble parameter H (right panel), as a function 
of redshift. The uncertainty is relative to a cosmological constant {wq = —1). 



MWA5000. The small change in the radial BAO scale with 
Wq is reflected by the large value of the derivative in the 
lower right panel of Figure [T] 



5.1 Constraints on parameterized models of 
evolving dark energy 

In the remainder of this section we constrain two different 
parameterized dark energy models assuming observations 
with the MWA5000 at four different redshifts. Rather than 
consider the spherically averaged constraint on /ca, in this 
section we utilize independent constraints of the line-of-sight 
and transverse BAO scales (Glazebrook & Blake 2005). Fig- 
ure [H] illustrates that the MWA5000 will have a compara- 
ble sensitivity to measurements of kA along the line-of-sight 
(A:a,||) as compared to perpendicular to it (A;a,±). We there- 
fore assume the transverse and line-of-sight sensitivities to 
be AfcA,x/fcA,x ~ AfcA,|| /A;a,|| = \/2 x AfcA/^A, and con- 
strain dark energy parameters using the joint constraints on 
Da and H (i.e. we assume AA;a,±/A;a,± = ADa/Da and 
AfcA,||/fcA,|| = AH/H). Note that our relations between the 
transverse or line-of-sight sensitivities and the spherically 
averaged sensitivity to the BAO scale assume that the trans- 
verse and line-of-sight scales are independent, whereas Fig- 
ure |6] illustrates that there is some degeneracy. We therefore 
slightly underestimate the uncertainties on the transverse 
and line-of-sight BAO scales. 

The equation of state of dark energy may not be con- 
stant. The possible evolution is often parameterized using 
the following simple form 



where wi is the derivative of wq (z) with respect to the scale- 
factor a. We have investigated the joint constraints that the 
21cm BAO observations might place on the parameters wq 
and wi , and present the results in Figure [S] We construct 
the relative likelihoods for wq and wi given a true model 



with wq 



C{WQ,Wl) 



-1 and wi = 0, 

1 
2 



exp < -- 



SkA,i 

"'^A,l 



+ 



5k 



A, II 



^iLobs 
A, II 



(15) 



WQ (z) = WQ + wi{l — a) , 



(14) 



where SkA,_L{z) and SkA^\\{z) are functions of wq and wi. 
We show results in 4 redshift bins and have assumed 
SkfX = Skf l = y^Skf", with observed uncertainties of 
Skf" = AkA/kA = 0.01 at z = 1.5, Skf"" = 0.007 at 
z = 2.5, 5fcA ' = 0.02 at z = 3.5, and Skf'' = 0.007 at 
2 = 6.5. The latter two cases are examples of the levels 
of precision in AfcA/^A that may be achievable in lOOOhr 
integrations over each of three fields with the MWA5000, 
while the cases at z ~ 1.5 and z ~ 2.5 correspond to an 
instrument with the antenna number of the MWA5000, but 
with a higher frequency range and = Atne (again over 
three fields and 1000 hours). In each case contours of the 
likelihood are shown at 64%, 26% and 10% of the maxi- 
mum. The contours illustrate the degeneracy in wo and wi 
from a single measurement of the acoustic scale at high red- 
shift. This degeneracy arises if the cosmological constant is 
the true model, because the effect of dark energy is limited 
to the value observed for Da- However both wq and wi 
could be constrained at z < 2.5. We note that models with 
wq{z) < — 1 or wq{z) > 1 violate the weak energy condi- 
tion, and shade these regions grey in Figure [S] However, we 
include constraints on all parameters for generality. 

Constraints on the parameterization of the dark en- 



© 0000 RAS, MNRAS 000, 000-000 



14 Wyithe, Loeb & Geil 




Figure 8. The joint lilcelihoods for wq and w\ relative to a model witli a cosmological constant {wq = —1, w\ = 0). The contours are at 
64%, 26% and 10% of the maximum and were found from combining constraints on Da and H, assuming uncertainties AA:a ±/kA ± = 
AfcAj|/A:A,|| = X 0.01 (at z = 1.5, upper left panel), A/cA.x/fcA.x = ^fcA.y/'^A.II = X 0.007 (at z = 2.5, upper right panel), 
AkA^±/kA,± = AfcA,||/fcA,|| = V2 X 0.02 (at z = 3.5, lower left panel), and AkA,±/kA,± = AfcA,||/fcA,|| = V2 X 0.007 (at z = 6.5, lower 
right panel). At z ^ 3.5 these sensitivities correspond to three fields, each observed for lOOOhr with the MWA5000. At z ^ 2.5 these 
sensitivities correspond to three fields, each observed for lOOOhr with an array having 5000 antennae (each with 16 cross-dipoles) and 
= ^tilc- Models with wq{z) < — 1 or wq{z) > 1 violate the weak energy condition, and we shade these regions grey. 



ergy in equation (|14|) are most accessible to observations at 
z < 2.5. Indeed by construction equation (|14p does not lead 
to differences in the expansion history ai z > 3.5 where the 
BAO scale is accessible to the MWA5000. Nevertheless, the 
combination of observations at several redshifts at z > 3.5 
would provide tight limits on the value of either wq (< 5%) 
if it is constant, or wi (±0.2) if wq is known apriori. In con- 
structing these limits we have utilized external information 
through our assumption that the cosmological parameters 
Jim, fib and JIq are known, but have not incorporated ad- 
ditional sources of constraint on the dark energy equation 
of state at lower redshift, which could include observations 
of Type la supernovae (e.g. Riess et al. 2007; Zhao et al. 
2007), or the BAO scale from galaxy redshift surveys (Cole 
et al. 2005; Eisenstein et al. 2005). The constraints on the 
dark energy equation of state presented in Figure [5] could of 
course be improved through a joint analysis of all available 
constraints. For example, while the BAO scale measured 
by Eisenstein et al. (2005) does not by itself constrain the 
parameters of dark energy, Wood-Vasey et al. (2007) have 
combined this measurement with observations of Type la 
supernovae from the ESSENCE supernova survey, and find 
(assuming a flat Universe) that wq = — 1.05to i2 with a 
systematic uncertainty of ±0.11. 

Equation (|14[) imposes a certain type of evolution on the 
dark energy. In particular, for all values of wi the parameter- 
ization in equation (|14|l assumes that 50% of the evolution 
in wq{z) occurred at z < 1. Moreover, the constraints on 
model parameters are only meaningful should the model be 
a correct description of reality. Since 21cm observations of 
BAO will be a powerful probe of dark energy at high red- 



shifts, we require a parameterization that is more general 
than is provided by equation (|14|) . Current observations indi- 
cate that the evolution of the cosmic expansion is consistent 
to within ~ 10% with a pure cosmological constant at 2 < 1 
(Riess et al. 2007; Zhao et al. 2007), while redshifted 21cm 
observations might be the first working method of probing 
the dark energy at 2; > 3. An alternative parameterization 
to equation (|14|) for the possible evolution of the vacuum 
energy density [pii(2)] at high redshifts is to adopt a con- 
stant value {pv,o) in the redshift range 0^2^ 2miii which 
is constrained by other observational probes (such as Type 
la supernovae or galaxy redshift surveys), and to represent 
the unknown evolution at higher redshifts 2 > 2min by 

pv = Pv,0 [1 ± Cl(2 — 2min) + (^2(2 — 2min)^ ± ...] , (16) 

where we have kept only the two leading terms in the poly- 
nomial expansion of pv{z). We next require that the Uni- 
verse be flat, which yields the corresponding evolution for 
the Hubble parameter at 2 > 2min, 
n 2 

= r2,„^o(l ± 2)^ + ^^t.,o[l ± Cl(2 — 2min) 



H{z) 



Ho 



±C2(2 



0^ + . 



(17) 



where the subscript denotes present-day values. This ex- 
pression can then be substituted directly into equation (|10[l 
to obtain the angular diameter distance, /3a (2)- If the vac- 
uum energy density originates from a scalar field </> rolling 
down a potential V(<?!>), then measurements of the above ex- 
pansion parameters can be used to infer the shape of V{(j)) 
(see Linder 2007 and Turner & Huterer 2007 for recent re- 
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Figure 9. The joint likelihoods for C\ and C2 relative to a model with a cosmological constant (Ci = C2 = 0). The contours are at 
64%, 26% and 10% of the maximum and were found from combining constraints on Dx and H, assuming uncertainties Ak\ j_/k\ j_ = 
^kA,\\/^A,\\ = X 0.007 (at z = 2.5, upper left panel), AA;A,x/fcA,± = ^'=A,||/'i^A,|| = X 0.02 (at z = 3.5, upper right panel), 
AkA,±/kA,± = ^fcA,||/'=A,|| = \/2 X 0.007 (at 2 = 6.5, lower left panel), and AfcA,_L/fcA,x = AfcAj|/A;A,|| = X 0.04 (at z = 12, lower 
right panel). At z ^ 3.5 these sensitivities correspond to three fields, each observed for lOOOhr with the MWA5000. At z = 2.5 these 
sensitivities correspond to three fields, each observed for lOOOhr with an array having 5000 antennae (each with 16 cross-dipoles) and 
= Atiic. Models with pt, < are shaded grey. 



We have investigated the joint constraints that the 21cm 
BAO observations might place on the parameters Ci and 
C2, and present the results in Figure O We do not constrain 
pv to be positive, but shade the regions where pv < 0- For 
this figure we assume Zmin = 1 and that the dark energy 
at z < 2min IS a cosmological constant. We construct the 
relative likelihoods for Ci and C2 given a true model with 
Ci = C2 = 0. As before we show results in four redshift 
bins and have assumed observed uncertainties of Sk'^^j^ = 
Skf l = V25kf'' with Skf'' = AkA/kA = 0.007 at z = 2.5, 
Skf'' = AkA/kA = 0.02 at z = 3.5, AkA/kA = 0.007 at 
z = 6.5, and AkA/kA ~ 0.04 at z — 12. In each case con- 
tours of the likelihood are shown at 64%, 26% and 10% of 
the maximum. Since this model does not exclude a dark 
energy contribution to the expansion dynamics at high red- 
shift, its parameters are well constrained by high redshift 
observations of the BAO scale. The contours iUustrate the 
level of degeneracy in C'l and C'2 from a single measurement 
of the acoustic scale at high redshift. The combination of 
constraints at different redshifts would reduce this degener- 
acy. 



6 DISCUSSION 

In this paper we have shown that measurements of the power 
spectrum of fluctuations in the intensity of 21cm emission 
could provide a precise determination of the scale of the 
acoustic peak in the matter power spectrum. By using the 
acoustic scale to measure the angular diameter distance and 
Hubble parameter, it would therefore be possible to use red- 



shifted 21cm studies to constrain the dark energy in the un- 
explored redshift range ofl.5<z<15. 

The main advantage of 21cm observations is that in 
addition to studies in the redshift range where a cosmologi- 
cal constant is expected to dominate the expansion dynam- 
ics, measurement of the BAO scale would also be possible 
at much earlier cosmic epochs than are accessible by other 
methods. Indeed, observations of the redshifted 21cm power 
spectrum at 3.5 < 2 < 12 using a facility with ten times the 
collecting area of the MWA would make a detailed study 
of the presence of dark energy at high redshift that differs 
from that measured at 2 < 1. Moreover, 21cm observations 
will be comparably sensitive to the line-of-sight and trans- 
verse acoustic scales. In contrast, spectroscopic galaxy red- 
shift surveys are more sensitive to the transverse scale, while 
photometric redshift surveys will not be sensitive to the ra- 
dial BAO at aU (Glazebrook & Blake 2005). 

We note that a measurement of the power spectrum 
both along the observer's line-of-sight as well as perpendic- 
ular to it can also be used to constrain cosmological param- 
eters through the Alcock-Paczynski (1979) test in a 21cm 
data set (Barkana & Loeb 2005a; Nusser 2005; Barkana 
2005). However, the use of the BAO scale as a cosmic yard- 
stick allows the BAO to be a more powerful probe of dark en- 
ergy than the Alcock-Paczynski (1979) test. This is because 
the line-of-sight and transverse measurements of the BAO 
scale yield independent measurements of the expansion his- 
tory, whereas the Alcock-Paczynski (1979) test utilizes only 
their ratio. 

The redshift range accessible to an instrument like the 
MWA is specifled by the antennae design. A purpose built 
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instrument with sensitivity to a higher range of frequencies 
could probe the BAO scale at lower redshifts, enabling con- 
tinuous coverage from the reionization era to 2 < 1 where 
galaxy and supernovae surveys become effective. Redshifted 
21cm studies of acoustic oscillations could in principle also 
be extended to the redshift regime 15 < 2: < 200 in the 
more distant future (Barkana & Loeb 2005b) . However, such 
an extension would be particularly challenging because of 
the increase in the Galactic synchrotron foreground towards 
lower observed frequencies. 

We have shown that the first generation of low- 
frequency experiments (such as MWA or LOFAR) will be 
able to constrain the scale of the acoustic oscillations to 
within ~ 3% in a high redshift window just prior to overlap 
(see Fig. [5]). This sensitivity to the acoustic peak is compa- 
rable to the best current measurements from galaxy surveys 
(Eisenstein et al. 2005). In addition, we find that observa- 
tions using an array 10 times the size of the MWA, and 
with sensitivity in an appropriate frequency range would 
achieve measurements of the BAO scale with precisions that 
are below one percent in some particular redshift intervals. 
These windows include the reionization era, as well as the 
range 1.5 < 2 < 3.5, where the precision obtained would be 
comparable to ambitious high redshift spectroscopic surveys 
of ~ 10® galaxies, covering more than 1000 square degrees 
(Glazebrook & Blake 2005). Moreover, since observations of 
the BAO scale in 21cm fiuctuations can be extended down 
to 2 ~ 1.5 there will be an opportunity to demonstrate that 
the results match those from other well studied techniques. 
However, at 2 > 3.5, where the acoustic peak scale is not ac- 
cessible through standard techniques, the 21cm power spec- 
trum would provide unique constraints on the dark energy. 

Acknowledgments The authors thank Chris Blake, 
Matt McQuinn and Brian Schmidt for comments which have 
greatly contributed to this work. The research was sup- 
ported in part by the Australian Research Council (JSBW) 
and Harvard University grants (AL). PMC acknowledges the 
support of an Australian Postgraduate Award. 



Furlanetto, S. R., Zaldarriaga, M., & Hernquist, L. 2004, 
ApJ, 613, 16 

Furlanetto, S. R., Oh, S. P., & Briggs, F. H. 2006, Phys. 
Rep., 433, 181 

Geil, P.M., Wyithe, J.S.B., 2007. l arXiv:0708.3716V l 
Glazebrook, K., & Blake, C. 2005, ApJ, 631, 1 
Hu, W., & Haiman, Z. 2003, Phys. Rev. D., 68, 063004 
Linder, E. V. 2007, to appear in General Relativity and 
Gravitation, e-print'arXiv;0704.2064v2 [astro- ph] 
McQuinn, M., Zahn, O., Zaldarriaga, M., Hernquist, L., & 
Furlanetto, S. R. 2006, ApJ, 653, 815 

Mesinger, A., & Furlanetto, S. 2007, ArXiv e-prints, 704, 
larXiv:0704.0946l 

Morales, M. F. 2005, ApJ, 619, 678 

Nusser, A. 2005, MNRAS, 364, 743 

Peacock, J. A., & West, M. J. 1992, MNRAS, 259, 494 

Percival, W. J., Cole, S., Eisenstein, D. J., Nichol, R. C, 

Peacock, J. A., Pope, A. C, & Szalay, A. S. 2007, ArXiv 

e-prints, 705, arXiv:0705.3323 

Prochaska, J. X., Herbert-Fort, S., & Wolfe, A. M. 2005, 
ApJ, 635, 123 

Riess, A. G., et al. 2007, ApJ, 659, 98 

Seo, H.-J., & Eisenstein, D. J. 2003, ApJ, 598, 720 

Seo, H.-J., & Eisenstein, D. J. 2005, ApJ, 633, 575 

Seo, H.-J., & Eisenstein, D. J. 2007, ApJ, 665, 14 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Turner, M. S., & Huterer, D. 2007, ArXiv astrophysics e- 

prints, arXiv:astro-ph/0706.2186 

Wood-Vasey, W. M., et al. 2007, ApJ, 666, 694 

Wyithe, J. S. B., & Loeb, A. 2004, Nature, 432, 194 

Wyithe, S., & Loeb, A. 2007, ArXiv Astrophysics e-prints, 

arXiv:astro-ph/0708.3392 

Zhao, G.-B., Xia, J.-Q., Li, H., Tao, C, Virey, J.-M., Zhu, 
Z.-H., & Zhang, X. 2007, Physics Letters B, 648, 8 



REFERENCES 

Abdalla, F. B., & Rawlings, S. 2005, MNRAS, 360, 27 
Alcock, C, & Paczynski, B. 1979, Nature, 281, 358 
Angulo, R., Baugh, C. M., Frenk, C. S., & Lacey, C. G. 
2007, ArXiv Astrophysics e-prints, ar Xiv:astro-ph/0702543| 
Barkana, R. 2006, MNRAS, 372, 259 
Barkana, R., & Loeb, A. 2005, ApJL, 624, L65 
Barkana, R., & Loeb, A. 2005a, MNRAS, 363, L36 
Bharadwaj, S., & Ah, S. S. 2005b, MNRAS, 356, 1519 
Blake, C, & Glazebrook, K. 2003, ApJ, 594, 665 
Bowman, J. D., Morales, M. F., & Hewitt, J. N. 2006, ApJ, 
638, 20 

Bowman, J. D., Morales, M. F., & Hewitt, J. N. 2007, ApJ, 
661, 1 

Cole, S., et al. 2005, MNRAS, 362, 505 

Corasaniti, P.-S., Huterer, D., & Melchiorri, A. 2007, Phys. 

Rev. D, 75, 062001 

Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 
Eisenstein et al. 2005, ApJ, 633, 560 



© 0000 RAS, MNRAS 000, 000-000 



